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EXECUTIVE  SUMMARY 


OBJECTIVE 

The  objective  of  this  work  was  to  examine  the  theoretical  formulation  of  shell  membrane 
waves  to  determine  how  a  computationally  efficient,  yet  faithful,  model  could  be  realized. 

RESULTS 

Faithful  replication  of  the  behavior  may  be  obtained  from  a  relatively  simple  set  of 
expressions.  These  produce  tremendous  computational  savings. 

RECOMMENDATIONS 

The  proposed  model  should  be  implemented  as  soon  as  feasible. 
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INTRODUCTION 


Accurate  target  models  capable  of  predicting  the  echo  time  history  of  submarines  and 
surface  ships  are  essential  tools  for  antisubmarine  warfare  and  ship  vulnerability  studies.  A 
realistic  model  must  account  for  the  many  highlights  of  a  target  as  well  as  the  conditions  of 
^  insonification.  The  model  must  predict  the  proper  highlight  amplitude  and  phase  as  well  as 

temporal  and  spectral  behavior.  Thus,  physical  scattering  mechanisms  must  be  properly 
identified  and  incorporated  into  the  model. 

•  Over  the  past  30  years,  target  models  developed  at  Space  and  Naval  Warfare  (SPAWAR) 

Systems  Center,  San  Diego  (SSC  San  Diego)  have  proven  to  be  an  invaluable  asset  for 
developing  signal-processing  algorithms  and  engagement  simulations.  These  models  have 
been  exported  throughout  the  Navy  community.  Previous  formulations  of  the  SSC  San  Diego 
target  models,  dictated  by  past  requirements  and  test  scenarios,  have  neglected  elastic  wave 
effects.  This  restriction  simplified  the  numerical  implementation  (as  well  as  theoretical 
development)  considerably.  This  restriction  serves  well  for  high-frequency  and  even 
narrowband  mid-frequency  applications.  However,  for  broadband  processing  schemes,  more 
sophistication  is  required. 

Efforts  have  been  undertaken,  e.g.,  interacting  ribs  (Lengua,  1997).  This  report  examines 
the  modeling  of  shell  membrane  waves.  We  will  begin  with  general  background  information 
since  many  interested  parties  are  unfamiliar  with  the  subject.  Then  we  will  discuss  the 
modeling  details. 


BACKGROUND 


MEMBRANE  SHELL  THEORY 

We  will  take  the  approach  of  beginning  with  the  simplest  shell  theory  and  then  proceed 
with  more  complicated  formulations.  The  theory  discussed  in  this  section  is  known  as 
membrane  shell  theory.  We  will  follow  the  development  of  Graff  (1975,  pp.  259-262).  Only 
forces,  both  normal  and  shear,  acting  in  the  midsurface  of  the  shell  are  considered.  The 
transverse  shear  forces  and  the  bending  and  twisting  moments  are  assumed  negligible.  Thus, 
the  shell  behaves  as  a  curved  membrane. 


Consider  a  cylindrical  shell  of  radius  a ,  thickness  h ,  and  density  p  .  We  will  use 
cylindrical  coordinates  (r,6,z).  The  displacement  components  are  w ,  v,and  u  in  the 
radial,  tangential,  and  longitudinal  directions.  Let  us  examine  the  forces  (membrane  stresses 
N  and  applied  load  q )  acting  on  a  differential  element  of  the  shell.  The  equations  of  motion 
in  the  longitudinal,  tangential,  and  radial  directions  are,  respectively. 
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The  membrane  stresses  {N ^ ,Ng,N^, N^,  )  are  obtained  by  integrating  the  usual  stresses 
(2^  ^ ^ t  )  across  the  shell  thickness.  From  Hooke’s  law,  we  have 
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where  and  e,  are  the  axial  and  tangential  strains  of  the  middle  surface  of  the  shell 
element  and  7  =  is  the  shear  strain.  E  is  Young’s  modulus  and  o  is  Poisson’s 

ratio.  For  the  assumption  of  membrane-type  stresses  only,  (2^,2^,^^^,^^.^)  are  constant 
across  the  shell  thickness.  So  then 

..Eh  ,  X 

N.  = - -{e^+oe,} 


N. 


"  1-0 
Eh 


^“1-0^ 


{£,+0£^) 


Eh 


2(1  +  0-)'  ■ 


Let  us  now  consider  the  dynamics  of  the  deformation.  For  the  present  conditions,  these  are 
relatively  simple.  In  the  axial  direction 
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The  expression  for  the  shear  strain  may  be  obtained  by  considering  small  changes  in  angle  of 
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The  equations  of  motion  may  finally  be  written  as 
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where  c. 


is  the  thin-plate  speed. 


From  this  last  equation,  we  see  that  the  normal  displacement  of  the  shell,  w ,  is  coupled  to 
the  two  “in-plane”  displacements,  u  and  v .  For  flat  plates,  the  normal  displacement  is 
independent  of  the  in-plane  displacements. 


3 


DONNELL  SHELL  THEORY 


Formulation 


The  analysis  of  a  shell  including  bending  effects  on  the  deformation,  as  well  as  bending 
moments  and  shear  forces  in  the  equations  of  motion,  yields  considerably  more  complex 
expressions  than  those  in  the  preceding  section.  The  Donnell  formulation  (Kraus,  1967, 
p.  297)  includes  simplified  versions  of  these  effects.  The  approximations  are  related  to  the 
influence  of  transverse  shear  forces  on  tangential  motion  and  to  the  expressions  for  curvature 
and  twist.  The  results  are  equations  very  much  like  the  membrane  equations  of  motion.  There 


is  only  an  additional  term  of 


h 

on  the  left  side  of  the  radial  displacement  equation 


w 

TT 


a" 


a 


dz^ 


w  ^  d*w  1  aw 
+  2.  „  + 


f-1 

3u 

1  'I 

dz 

3v  1-(T^ 

- j - q  ; 

dd  Eh 


1  a^w 


where  = 


12a" 


We  wish  to  study  the  propagation  of  harmonic  waves  in  the  shell.  These  will  have  an 
angular  frequency  (a  and  wavenumber  k .  Note  that,  since  we  are  using  a  “thin-shell” 
analysis,  it  is  required  that  kh«n  ,  or  equivalently  pka  «  1 .  For  high  frequencies,  a  thick- 
plate  analysis  will  be  necessary  (Junger  and  Feit,  1986,  pp.  214—215).  Rather  than  begin  with 
the  general  solution,  it  is  illustrative  to  consider  some  special  cases. 


Case  of  0-Independent  Motion 

One  of  the  most  important  special  cases  arises  from  considering  motion  independent  of  6  . 
Here  the  equations  of  motion  are  (also  set  q  —  Q) 
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Note  that  the  equation  for  tangential  motion  has  uncoupled  from  the  remaining  equations. 
It  may  be  written  as 
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where  = 


is  the  shear  speed.  This  is  a  one-dimensional  equation  describing  the  purely 


torsional  motion  of  the  shell.  The  propagation  speed  is  the  same  as  for  such  waves  in  a  solid 
circular  rod  (Graff,  1975,  p.  263). 


Consider  now  the  coupled  equations  in  u  and  w .  Let  u  -  and  w  =  A 
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The  determinant  of  coefficients  gives  the  frequency  equation 


where  Q  =  —  ty  is  the  normalized  frequency.  In  terms  of  the  phase  velocity  c  =  — ,  we  have 
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The  long-  and  short-wavelength  limits  are  easily  obtained.  At  long  wavelengths  (ka  -^0) 
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so  one  solution  is  c  ^  j —  =  ,  the  longitudinal  bar  speed  (Graff,  1975,  p.  264).  In  this 
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W 


u 


c  — .  Here 
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ka,  so  the  motion  is  primarily  longitudinal  (as  expected).  Another  solution  is 
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At  short  wavelengths  {ka-^o°) 
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so  one  solution  is  c  ->  .  Then  F  =  0 ,  the  motion  is  primarily  longitudinal  (quasi- 


U_ 

W 


1 


oc  —  .the  motion  is 
ka 


compressional  wave).  Another  solution  ffcaCp .  Here 

primarily  radial  (quasi-flexural  wave). 

Figure  1  shows  the  dispersion  curves  when  cr  =  0.3  and  =  0.01 .  Figure  2  shows  |t/| 
fijr  the  two  branches  (normalized  such  that  |t/j^  +  \wf  =1). 


Figure  1 .  Example  of  dispersion  curves  for  0-independent  motion. 
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Figure  2.  Example  of  longitudinal  displacement  for  0-independent  motion. 


Case  of  z-Independent  Motion 

Another  important  special  case  arises  from  considering  motion  independent  of  z .  Here 
the  equations  of  motion  are  (also  set  =  0 ) 
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Note  that  the  equation  for  longitudinal  motion  has  uncoupled  from  the  remaining  equations. 
It  may  be  written  as 

b^u  ^  1  b^u 

b{a9Y 

so  that  again  disturbances  propagate  at  the  shear  speed. 

Consider  now  the  coupled  equations  in  v  and  w .  Let  v  =  and  w  =  . 

Now  we  must  have  our  first  discussion  of  boundary  conditions.  These  are  namely  continuity 
of  the  displacements  for  9  and  9  + In 

v(^  +  2;r)=  v(^) 
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and  similarly  for  w .  Thus  m-n,  where  n  is  an  integer  (>  0).  Therefore,  only  discrete 
modes  exist.  Note  that  we  must  have  pn  « 1  for  the  development  to  be  valid. 


In  terms  of  the  phase  velocity 

2 

For  the  «  =  0  mode,  one  solution  is  c  ^  0 .  Here  co  —  0  and  W  =  0 ,  so  it  is  a  trivial 

c 

solution.  Another  solution  is  c  — >  o® .  In  this  case,  co  =  —  and  V  =  0,  which  is  a  pulsating 

a 

cylinder. 

At  short  wavelengths  ( n  — >  o© ),  one  solution  is  c  ^  .  Then  W  =  0 ,  the  motion  is 

primarily  tangential  (quasi-compressional  wave).  Another  solution  is  c  — >  pnc^ .  Here, 

oc  — ,  the  motion  is  primarily  radial  (quasi-flexural  wave). 
n 

Figure  3  shows  the  dispersion  curves  when  p  =  0.01 .  These  are  plotted  as  continuous 
functions  of  n  as  a  guide.  Figure  4  shows  |f|  for  the  two  branches  (normalized  such  that 

!vr+M^=i). 
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Figure  3.  Example  of  dispersion  curves  for  z-independent  motion. 


Figure  4.  Example  of  tangential  displacement  for  z-independent  motion. 


Forced  Planar  Vibrations 

Let  us  again  consider  z-independent  motion  but  now  with  an  applied  load  that  is 
symmetric  in  9 . 
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n=l 

The  equations  of  motion,  after  using  the  trigonometric  orthogonality  relations 
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These  are  divergent  at  the  modal  frequencies  because  dissipation  has  been  ignored. 

Later  it  will  be  convenient  to  express  the  response  in  terms  of  the  modal  mechanical 
impedance 
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so  that  is  the  ratio  of  the  modal  pressure  on  the  outer  shell  surface  to  the  modal  radial 
velocity. 

Note  that  we  have  been  considering  the  in  vacuo  vibrations  of  the  shell.  A  shell  vibrating 
in  a  fluid  will  radiate  into  the  fluid.  The  radiation  loading  is  significant  and  will  be  discussed 
later. 


General  Motion 

In  considering  general  motion,  let 

V  = 

Solution  of  the  equations  of  motion  (with  q  =  0)  requires 
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The  determinant  of  coefficients  gives  the  frequency  equation 


where  the  parameters  are 
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and  we  have  defined  a  helical  wavenumber 


c  Q 

The  phase  velocity  is  then  —  = - .  Here,  we  require  pk^a  «  1  for  the  development  to  be 

Cp  A" 

valid. 

For  the  n  =  0  mode,  we  retrieve  the  0-independent  results.  Similarly,  if  A:  =  0 ,  we 
retrieve  the  z-independent  results. 

Let  us  consider  the  situation  when  A:^a  »  1 .  The  phase  velocity  equation  may  be  written 
as 
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The  solutions  are  then  c  -»  ,  c  ^  pks^^p ,  and  c  ^  c, .  To  examine  the  displacements, 

we  must  know  the  relative  contributions  of  ka  and  n .  Let  us  first  suppose  ka»  n  .  For 
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,  the  motion  is 
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(quasi-flexural  wave).  For  c  c^, 

tangential.  Now  suppose  n»  ka  .  For  c  c^, 

primarily  tangential  (quasi-compressional  wave).  For  c  , 

and  the  motion  is  again  primarily  radial  (quasi-flexural  wave).  For  c  , 
W  ka 

—  oc  — ,  the  motion  is  primarily  longitudinal. 

U  n 


The  frequency  equation  is  a  cubic  equation  in  and  may  therefore  be  solved  using  the 
standard  method  (Gautschi,  1965,  p.  17)].  Let 

1  .  1  .2 


u 

ka  , 

y 

oc  — -  and 

w 

w 

u 


n 


ka 

—  and 


<2  =  -A — K 

3  9 

«  =  -|(A4-3A,)+^A’  ■ 


If  +R^  >  0 ,  there  is  one  real  root  and  a  complex  conjugate  pair.  H  +R^  =0,2M 
roots  are  real  and  at  least  two  are  equal.  If  +R^  <  0 ,  all  roots  are  real  and  distinct.  Let 


*^1,2  ~ 


1/ 


Then 


nf=s,+s,+i/i, 


£2?  =-i(5, +S,)+i4  +^i{S, -s,) 


Once  the  frequencies  have  been  evaluated,  the  phase  velocities  can  easily  be  determined. 
The  displacements  may  be  calculated  from  the  expressions 
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U^aha 

2  2(7 

a^-(fezy-i(i-<Ty  ^ 

—  =:  i - ± - i - n — 

U  oka  2(T  U 

along  with  the  normalization  jC/[^  +  jF|  +  |fF|  =  1 . 

Let  us  consider  some  examples.  We  will  again  let  rr  =  0.3  and  ^  =  0.0 1 .  It  is  useM  to 
review  the  n  =  0  case.  Figure  5  shows  the  frequency  roots,  while  figure  6  shows  the 
corresponding  phase  velocities.  Note  the  crossover  of  roots  2  and  3.  The  displacement 
amplitudes  for  the  three  roots  are  shown  in  figures  7  through  9.  The  tangential  displacement 
does  not  appear  because  it  is  decoupled  from  the  others.  The  behavior  of  the  n  =  1  case  is 
shown  in  figures  10  through  14,  while  that  of  the  «  =  2  case  is  shown  in  figures  15 
through  19. 
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Figure  6.  Phase  velocities  for  the  n  =  0  case. 


Figure  7.  Displacements  for  root  1  of  the  n  =  0  case. 


MAQNITUDH 


0  0.5  1  1.5  2  2.5  3  3.5  4  4.5  5 

ka 

Figure  9.  Displacements  for  root  3  of  the  n  =  0  case. 
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Figure  1 1 .  Phase  velocities  for  the  n  =  1  case. 


Figure  14.  Displacements  for  root  3  of  the  n  =  1  case 


Figure  1 5.  Frequencies  for  the  n  =  2  case, 


Forced  Vibrations  When  |3  =  0 

Let  us  consider  forced  vibrations  from  an  applied  load,  which  again  is  symmetric  in  6  but 
now  has  a  harmonic  variation  in  z 

q  = 

with 

/(^)=X/nCOs(rt^). 

n=0 

For  simplicity,  we  will  assume  fi  =  0. 

Let 

u  =  ^U„  cos(n^>'^*^-‘“^ 

«=0 
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Here, 


A  = 


1-0 


+  (k^af  +  2(1  +  aXkaf  ] 


It  may  be  easily  verified  that,  as  k  —>0 ,  u  —>0  and  our  earlier  results  for  v  and  w  are 
recovered. 


Thus,  the  modal  mechanical  impedance  is 

pCph 


z:{n,kah-i- 


Q  a 


and 


3m  .  , 

—  =  -laka  > 


2(7  ""  /.cos(ne>'<*'-“) 


n=l 


b’-CM/] 


COs(«g>'^^ 

3t  z! 


Again,  note  that  these  are  the  in  vacuo  vibrations  of  the  shell. 
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FINITE-LENGTH  SHELLS 


For  a  cylindrical  shell  of  finite  length,  the  specification  of  boundary  conditions  is  needed.  The 
simplest  set  is  that  of  a  simply  supported  shell.  If  the  ends  are  located  at  z  =  ±— L ,  the  displacement 

boundary  conditions  are  (lunger  and  Feit,  1986,  p.  218) 

3m  1 

w  =  — ^  =  v  =  — =  0  for  z  =  ±— L. 
dz^  dz  2 

The  motion  of  the  shell  may  be  described  by 

«  =  cos(n^)sin(/r„z>'““ 


V  =  sin(n^)cos(A:^z>-'‘“ 

m,n 

cos(n^)cos(^„z>‘"“ 

m,n 


where  =  (2m  + 1)— .  A  solution  with  the  cos(n6  )  and  sin(n6  )  factors  interchanged  would  also 
L 

be  valid. 

For  reference,  the  boundary  conditions  for  a  shell  that  is  clamped  at  its  ends  are 

u  —  v  =  w  =  ^^  =  0  for  z  =  ±— L. 

dz  2 


RADIATION  LOADING 

A  submerged  shell  undergoing  vibrations  will  radiate  into  the  surrounding  fluid.  As  usual,  we 
will  assume  a  time  dependence  e~'“‘ .  Let  P(r,6,z)  denote  the  radiated  pressure  field,  and  let  be 
the  density  of  the  fluid  and  Cg  be  the  wave  speed  in  the  fluid.  In  the  fluid,  P  must  satisfy  the 
Helmholtz  equation 


where  k  =  —  .In  cylindrical  coordinates,  this  is 


3^  1 

dr^  r  dr  dd^ 


+ 


P^k^P 


0. 


The  standard  solution  is  to  use  separation  of  variables.  Let  P(r,6,z)=  P(?')0(6)z(z).  Then 
Z  =  0  =  e""^and 
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dr  r  dr  r 


R  =  0. 


Therefore,  R  =  A7„  {jk^  -  r)+  fiF„  {^k^  -  k^  r),  where  7„  (x)  and  F„  (x)  are  Bessel  and  Weber 

functions  (Gautschi,  1965,  p.  358). 

The  boundary  condition  relating  P  to  the  shell  motion  is 


dr 


-Po 


This  represents  continuity  of  the  normal  component  of  the  force  at  the  boundary.  Since  a  fluid  does 
not  support  shear  motion,  the  tangential  components  of  the  force  may  be  discontinuous. 

The  radial  displacement  will  be  written  as  in  the  previous  section.  Since  we  desire  a  solution  in 
the  form  of  an  outgoing  wave,  the  pressure  field  is  expressible  as 

P{r,e,z,t)=  r)cos{n9)cos{k^zy““ 

m,n 

where  is  the  Hankel  function  of  the  first  kind  (Gautschi,  1965,  p.  358).  is  determined 

from  the  radial  boundary  condition 


mn 


_ PqW^ _ 


where  =  -CO^W^^ ,  and  the  prime  denotes  differentiation  with  respect  to  the  argument. 

Note  that  the  surface  pressure  obtained  from  these  equations  may  be  written  as 

Fia.e,z,t)=  cos(n6l)cos(*:.z>-" 


where 


Z'  = 

mn 


iPpC^k 


is  the  modal  radiation  impedance.  Thus,  no  energy  is  radiated  when  k  <k^. 

The  Hankel  functions  may  be  approximated  by  their  Debye  asymptotic  expansions  (Gautschi, 
1965,  p.  366). 


I  sin  7 


1  +  i 


,  5  2 

l  +  -cot  r 
3 

8z  sin  Y 


±i  z(siny-ycosr)-— 

,  L  4j 
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V  nz 


l  +  3i 


1  7  2 

1  +— cot  y 

9 

8zsin/ 


±i  z(sin7-f COS/)-- 

7 


in  which  y  =  cos  *f— 1  and  0  <  Re(7  )<  n  .  These  expansions  are  generally  valid  for  |z|  >  |v|  and 


\}i 


z  —  y  >  V  .  Note  that  sin  y  = 


.  Thus, 


- 1 . 
+  — z- 


It  is  worth  commenting  that  some  authors  ignore  the  cot^  y  term  in  the  Debye  expansions  and 
consequently  have  an  invalid  approximation  for  the  ratio  (Rumerman,  1996). 


Now,  if  »  n , 


^mn  ~  Po^O  ■ 


The  case  of  z-independent  motion  may  be  obtained  by  simply  setting  =  0,  Then, 

.  HP(ka) 

~  ^Po^O  TTTfi 


H’»(ka) 


and,  if  fai  »  n , 


the  specific  acoustic  impedance  of  the  fluid. 


Po^O 


The  effect  of  radiation  loading  on  forced  vibrations  may  be  easily  determined.  From  the  equation 
for  in  vacuo  vibrations  (generalized  to  axial  wavenumbers  .  In  a  fluid,  /„„  is 

replaced  by  ,  where  the  minus  sign  accounts  for  the  direction  of  the  force.  Now,  as 

discussed  earlier,  =  Z^„W„„  so  that 

’  mn  mn  mn 


W^n  = 

mn 


f  -Z’^  w 

J  mn  mn  mn 


or 


W  = 

mn 


fm 


2“  +Z'' 

mn  mn 
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Thus,  the  total  (or  shell-fluid  system)  modal  impedance  is  the  sum  of  the  shell  and  radiation 
modal  impedances,  as  might  have  been  expected. 


SCATTERING  OF  PLANE  WAVES 
Normal  Incidence 

Incident  Wave.  Let  us  consider  a  plane  acoustic  wave,  given  by  P^  {x,  y,  t)  =  Pq  exp[j (fex  -  tot)] , 
normally  incident  upon  an  infinite  cylindrical  shell  in  a  fluid.  In  cylindrical  coordinates 

P,{r,d,t)= 

Note  that  6=0  corresponds  to  the  forward-scattering  direction.  We  may  represent  P^  as  a  Fourier 
expansion  (Gautschi,  1965,  p.  361). 

P,  {r,e,t)=  Poe~““'^e„i"J„  {kr)cos{n0). 

n=0 

Note  that  since  (z)=  (- 1)"  (z),  we  may  write 


Pi{r,0,t)=  PqC  J ^{kr)cos{n0) 

n=-oo 


and  eliminate  the  factor.  It  is  often  useful  to  use  a  representation  in  terms  of  Hankel  functions 

cos(ne)[Hf(*r)+ 

^  n=—oo 


Rigid  Body:  Blocked  Pressure.  In  analyses  of  scattering  from  elastic  shells,  it  is  useful  to  first 
consider  scattering  from  a  rigid  object  of  the  same  shape.  The  scattered  pressure  field  will  be 
denoted  .  The  total  pressure,  which  in  this  case  is  termed  the  blocked  pressure,  is  then 

Pt  =P+P“®'‘*- 

D  I  S 

Since  the  body  is  rigid,  the  resultant  fluid  acceleration  must  have  a  zero  component  along  the 
normal  to  the  boundary 


at^ 


av, 


at^ 


=  0. 


a^ 


at^ 


is  the  normal  fluid  acceleration  that  would  be  observed  on  the  surface  in  the  absence  of  the 


scatterer,  that  is,  when  the  pressure  field  is  identical  to  the  incident  pressure.  This  is  given  by 


Po' 


av, 


at^ 


dp 


dr 


Combining  these  equations  gives 
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2„, Rigid 


Po- 


de 


dP. 


dr 


The  corresponding  surface  acceleration  distribution  is 


-w 


(6>)=  Po  — (*a)cos(nP). 


Po  n=0 


So, 


W„=P,—£XJ:{ka). 

Po 

Substituting  these  coefficients  into  the  expressions  from  the  Radiation  Loading  section  gives 

2'(fai) 


P,‘^{r,e,th-P„e-“'Ze/ 


n^O 


H'}'\ka) 


■Hf^{fcr^cos{n6). 


Thus,  the  blocked  pressure  is 

n  =  Pce'-i.e.r 

n-0 

or  alternatively 

n  “In*"”  i  >'■ 


cos 


(ne) 


cos 


(ne). 


The  subscript  1  corresponds  to  outgoing  waves,  and  the  subscript  2  corresponds  to  incoming  waves. 

Elastic  Shell.  While  a  rigid  scatterer  distorts  the  pressure  field  by  interfering  with  the  propagation  of 
the  incident  wave,  the  dynamic  response  of  an  elastic  scatterer  further  modifies  the  pressure  field. 

The  pressure  scattered  by  an  elastic  body  will  be  denoted  .  The  total  pressure  is  then 

Pt  =Pi 

It  is  convenient  to  express  Pf^'  as  the  sum  of  Pf^^^  and  an  unknown  component  P^  (lunger  and 
Feit,  1986,  p.  343) 

pElast  _  p  Rigid  p 


To  interpret  ,  note  that  the  boundary  condition  (on  radial  motion)  is  satisfied  if 


3P, 


dr 

dp 


dr 


d^w 

0p  Rigid 


dr 
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The  latter  condition  is  automatically  fulfilled  through  the  definition  of  .  The  former  condition 
indicates  that  is  equivalent  to  the  radiated  pressure  field  due  to  the  acceleration  of  the  elastic  body 
responding  to  the  blocked  pressure  field. 

We  will  use  the  Hankel  function  representation  of  the  blocked  pressure  to  allow  comparison  with 
the  results  of  Rumerman  (1991).  We  therefore  seek  a  solution  in  the  form 

PXr,0,t)=  ^  (kr  )cos(n  0^'"" . 


The  are  determined  from  the  boundary  condition 

PoW„ 


A„  = 


or  equivalently 


A„  = 


kH'J»{ka) 


iPo^^O^n 


Now 


bn 


w„  =  — 

zi+z: 


where  the  minus  sign  accounts  for  the  direction  of  the  force,  and  again  and  Z'  are  the  modal 
shell  and  radiation  impedances.  Consequently, 


A  =--P 


2  ''z:+z:H'}\ka) 

Thus,  the  total  pressure  field  may  be  represented  by  the  normal  mode  series 


Xi"  cos(ne)  H^\kr)- Hi\kr)R.(aM'>-^ 


HT\ka) 


Hy(ka) 


where 


R„{fl,ka)=' 


Zn  zt'{1)(u  \ 

_ Hn  \ka) 

r„{Q,ka)+iPoCo 

Hn\ka) 


/?„  is  a  generalization  of  the  plane- wave  reflection  coefficient,  to  which  it  reduces  when  ka  »  1  and 
ka»  n . 


The  denominator  of  /?„  is  the  modal  impedance  of  the  shell-fluid  system 
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Zr{Cl,ka)=  Z'^{Cl,ka)+  iPi,c„ 


Note  that  the  shell  modal  impedance  may  be  written  as 


Z‘(Q.,ka)=-ipc—Q. 


1  — 


Here,  Q.-~ka . 
c„ 


Rumerman  (1991)  expresses  the  total  pressure  field  as  a  contour  integral  with  n  generalized  to  a 
complex  number  v  .  This  is  noteworthy  because  the  poles  of  the  integrand  are  the  zeroes  of  the 
system  impedance,  viewed  as  a  function  of  v  .  Since  is  an  even  function  of  v  ,  the  poles  appear 
in  equal  and  opposite  pairs.  There  is  a  quasi-real  pole  with  Re(v  )=  Q  corresponding  to  a  quasi- 

compressional  wave.  Another  quasi-real  pole  with  Re(v')=  ■\l^l  corresponds  to  a  quasi-flexural 
wave.  These  are  generalizations  of  our  previous  results,  modified  by  the  effect  of  radiation  loading. 
However,  here  there  is  an  infinity  of  poles  v  creep  corresponding  to  creeping  waves.  There  are  also 
poles  corresponding  to  the  exponentially  decaying  quasi-flexural  near-field. 


Rumerman  shows  that,  when  |Vpj^^  \«ka,  the  flexural  wave  is  very  poorly  coupled  to  the  fluid 

and  makes  an  insignificant  contribution  to  the  scattered  field.  The  creeping  waves,  though  important 
in  forward  scattering,  decay  so  rapidly  as  they  circumnavigate  the  shell  that  their  backscattering 

contributions  are  also  small.  Thus,  if  we  restrict  ourselves  to  the  backscattering  half-space  |0|  >  , 

these  terms  may  be  neglected.  An  additional  simplification  is  that  quite  often  « 1 ,  and  when 
|v|  <  fcz  ,  quantities  relating  to  may  be  accurately  determined  with  the  bending  term  ignored. 
Therefore,  we  may  take 


h 


Zl  (Q,  fai)  «  -ipc  —  Q 


a 


1  — 


1 


Now,  in  the  backscattering  half-space,  for  most  or  all  of  that  space  (depending  on  the  value  of  ka  ), 
the  Hankel  functions  may  be  approximated  by  their  Debye  asymptotic  expansions.  Then, 


Z:{Q,ka)- 

Wo 

0 

0 

1 

2kasm‘^  y 

sin/ 

zya,ka)- 

,  Po^o 

2kas,m^  y 

~Y 

sin  7 

and 


h 


Zl(Q.,ka)^-ipCp-Q 


a 


■n 


■  Po^o  I  Po^o 
2kasm'^  y  sin 7 
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^  .  ^{kaf-v^ 

where  sin  y  = - 

ka 


The  system  pole  V^onip  is  found  by  setting  -  0.  An  approximate,  though  very  accurate, 
solution  is  found  by  letting  sin  j  =1 


'^Comp  ~  ^ 


-1-1 


1  +  _L£2." 

ka  p  h 


1  — 


2ka 


yj 


We,  of  course,  are  interested  in  values  of  ka  where  is  quasi-real  and  an  integer.  That  is,  we 
want  Vcomp  =  n(l  +  )  with  |^|  « 1 .  Here, 


(  \ 


) 


w- 


2  p  h 


{kaf  + 


p^a 

p  h 


V 


1  ka 


p^a 

p  h 


2n^ 


{kaf  + 


f  \2 

p  h 


Note  that  these  conditions  cannot  be  satisfied  by  n  =  0 .  The  n-0  mode  does  not  couple  to  the  fluid 
and  may  be  ignored  (Rumerman,  1993).  More  fundamentally,  an  imaginary  part  of  the  root  is  needed 
to  couple  to  the  fluid  (Rumerman,  1992).  For  typical  values  of  the  parameters,  the  condition  |(^|  « 1 

is  satisfied  if  ita  >  3/2 .  A  check  shows  that  siny-^l-(co/Cpf  ~1,  as  was  earlier  assumed. 


Oblique  Incidence 

Let  us  now  consider  a  plane  acoustic  wave  obliquely  incident  upon  an  infinite  cylindrical  shell  in 
a  fluid 

Pi(r,6,z,t)=  Pq  exp[/ (kr cos 6  cos ^  -I-  kz sin cut)] 

where  ^  is  the  angle  of  incidence  (with  respect  to  the  normal).  This  problem  is  formally  identical  to 
that  of  normal  incidence  through  the  replacement  of  k  by  K  =  k  cos  (j)  (or  of  the  sound  speed  in  the 
fluid  Cq  by  Cq/cos^  ).  However,  the  resulting  behavior  has  significant  differences  (Rumerman, 

1992).  These  are  that  an  obliquely  incident  wave  may  excite  two  kinds  of  membrane  modes  on  the 
shell,  while  only  one  is  excited  at  normal  incidence,  and  that  each  of  these  modes  is  a  supersonic 
wave  only  within  a  range  of  angles  about  normal  (and  evanescent  elsewhere). 

The  incident  pressure  may  be  expressed  as 


31 


P.  {r,e,  Z,r)=  cos{netH^\Kr)+  Hf\Kri 


so  that  the  blocked  pressure  is 

where  we  have  defined  an  axial  wavenumber  k^=k  sin  <p  . 

In  a  manner  analogous  to  that  used  in  the  previous  section,  the  total  pressure  may  be  represented 
by  the  normal  mode  series 


Jr  cos(npj 

^  ff=— <90 


where 


R^{Q.,ka,^)  =  - 


^  cosip  H'J-^\Ka) 
^  cos4>  H'„^^\Ka) 


The  modal  system  impedance  is 


Zj.  (Q,  ka,  k^a)=  Z*  (Q,  k^a)+  i 


o,c,  Hi\Ka) 


cos^  H'}\Ka)' 


The  total  pressure  may  be  expressed  as  a  contour  integral  (Rumerman,  1993).  As  before,  the  poles 
appear  in  equal  and  opposite  pairs.  Here  too,  when  \«ka,  the  flexural  wave  is  very  poorly 

coupled  to  the  fluid.  If  we  again  restrict  ourselves  to  the  backscattering  half-space  |P|  >  /r/2 ,  the 
creeping  wave  contribution  may  be  neglected. 

We  will  again  ignore  the  bending  terms  in  the  modal  shell  impedance  (Rumerman,  1993)  which 
may  be  written  as 

^R{k„ay  +Bj,{k^af  H-C^ 


z;(ak„a)=-/ 


pc„  h 


L  ^  ^\DAk.ciy+E,{k^af  +  F, 


with 


1-0 

2 

1-0 


[q'-(i-o-")] 


P,  =  (1  -  ~ - 1]-  [q^  -  (l  -  o-^ )] 
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Here  we  will  not  be  able  to  make  the  simplifying  approximation  sin  J  =1  to  find  the  poles. 
Instead,  replace  V  with  Ka cos}  ,  then  replace  cos^  7  with  1  -  sin ^  7,  and  solve  for  sinj  .  This 
yields  six  independent  solutions,  most  of  which  are  discarded  because  they  violate  the  Debye 
expansion  criteria.  The  two  valid  solutions  are  denoted  y  and  }  51,^3^ .  The  corresponding  poles 

are  Vco„p  =  .fiTacos^comp  and  =  Kacosy^i,,^ .  Rumerman  (1993)  discusses  the  typical 
behavior  of  these  poles  for  kacos^  >  3/2 .  When  (/>  <  ^^omp  =  sin"'  (cg/c^  ),  v^omp  is  essentially  real 
(with  a  small  imaginary  part)  representing  a  supersonic  wave.  When  ^  ,  Vcomp  is  essentially 

imaginary  (with  a  small  real  part)  representing  an  evanescent  mode.  Similarly,  when 
(p  <  =  sin  (cq  /Cj ) ,  V  is  essentially  real  representing  a  supersonic  wave,  and  when 

^  ^  ^shear  >  ^  Shear  is  esscntially  imaginary  representing  an  evanescent  mode. 

Since 

sin  7  = 
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COS0  = 


the  modal  radiation  impedance  may  be  expressed  as 

1 


Zl{ka,k^a)^  p^c^ka 


--1 


l!-MV 


Thus,  the  zeroes  of  the  modal  system  impedance  may  be  found  from 


0  =  -/ 


P^h 

Q  a 


-klV-y^f  +EAk.af 


Note  that,  if  =  0 ,  this  becomes 

'  pCp  h 


0  =  — / 


\kaf 


2  c 

- Q.  —  —ka  is  purely  real.  The  associated  membrane  mode  has  purely  in-plane 

'l-<7  c. 


Then,  Vs^ear  =  ■ 

motion  that  does  not  couple  to  the  fluid  (Rumerman,  1992).  It  may,  therefore,  be  disregarded.  Now, 


1- 


^_V  V 

ko-  j 


>=  1 ,  and  the  previous  result  for  v^omp  is  retrieved. 


The  equation  for  the  zeroes  of  the  modal  system  impedance  must  be  solved  numerically,  and  is 
tedious.  Figure  20  shows  the  first  five  compressional  and  shear  modes  when  a  =  3  m,  h  =  0.05  m, 
/O  =7800kg/m^,  CJ  =0.3,  £  =  2.04x10“ N/m^  p^  =1000  kg/m^  and  Cq  =1500 m/s.  Here, 

Cp  =5360m/sand  =3170m/s,  so  ^co„p  =16.3°  and  =  28.2°. 
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Figure  20.  Example  of  dispersion  curves  as  functions  of  frequency 
and  incident  angle. 


AMPLITUDE 

Rumerman  (1993)  has  shown  that,  by  use  of  a  Sommerfeld-Watson  type  transformation,  the 
scattered  field  may  be  expressed  in  terms  of  a  series  of  axial  modes  having  coefficients  that  are 
waveforms  in  the  circumferential  direction.  The  remainder  of  the  field  is  taken  as  a  “geometric” 
contribution. 

F{R,0J)=P^  +F" 

where  are  the  spherical  coordinates  of  the  observation  point. 


We  are  interested  in  the  wave  component  of  the  scattered  field.  In  the  far  field,  this  may  be 
evaluated  using  the  method  of  stationary  phase.  For  simplicity,  the  details  will  be  omitted.  The 
backscatter  result  is 


<s) 
■  m 


where  the  parenthetical  superscripts  distinguish  between  the  compressional  and  shear  contributions, 
and  is  the  largest  integer  not  exceeding  (l  +  while  is  the  largest  integer  not 

exceeding  -¥0X110^)117: .  Here,  the  normahzed  length  L  =  Lja .  Recall  that  K„^^a  =  . 

The  parameters  associated  with  the  incident  field  have  been  explicitly  called  out  for  clarity.  The 
coefficients  are  given  by 
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T  C  _  C  T 

^ia^L-  ^  ^ia„L 

1  _ 

1 

T  C  _  C  T 

_^P„v‘^a„v  ‘^p^v^a^v  _ 

,Zv(aa™)_ 

_  mK-a-^I  _ 

sin(v’;r) — ^ 
dv 

^-2iK^a(sia  S„  -S„  cosS„ ) 


sin^„ 

m 


Here,  and  are  the  two  linearly  independent  nondimensional  axial  wavenumbers  in  the 
membrane  range  (for  brevity,  the  superscripts  j  have  been  omitted),  =  ka  ,  and 
cos<5„,  ^v/{K  incfl).  The  parameters  S  and  T  are  normalized  displacements  and  stresses.  Their 
detailed  expressions  are  not  needed  in  what  is  to  follow. 

A  particular  mode  will  contribute  strongly  to  backscattering  if  two  conditions  are  satisfied.  The 
first  is  that  Re(v„  )  is  approximately  equal  to  an  integer  and  Im(v„, )« 1 .  This  condition  is 

independent  of  the  angle  and  signifies  the  mode  is  close  to  resonance  at  the  frequency 

considered.  The  second  condition  is  that  =  fajsin^i^  is  approximately  equal  to  a^.  This 

indicates  the  axial  variation  of  the  projection  of  the  incident  plane  wave  on  the  shell  is  close  to  the 
axial  variation  of  the  mode  shape. 

Because  Zj’ )  =  0  for  v  =  v„, ,  (a^  - a„  YJ )  when  , 

where  the  prime  denotes  differentiation  with  respect  to  a^.  In  addition, 

Z'J  Xda„/dv)=  -Zj  (iQ,a„  )  where  the  dot  denotes  differentiation  with  respect  to  v  . 
Therefore  when  ~ 


x  C  _  C  T" 

‘^p„v^a-^v 

(-ifi.  1 

_  "^P^v^a^v  ~  ^P„v1'a„v  _ 

_  mTC-aijl  _ 

sin(v;r)zJ(a,or„,) 

^-2iK-^ais.m  S„  -S„  cos  S„  ) 


L  sin^„  J 

The  quantity  in  the  second  set  of  brackets  is  a  directivity  function  of  that  mode.  If  the  angular 
width  of  the  mode  is  less  than  the  angular  distance  to  the  adjacent  modes,  the  peak  level  of  a  resonant 
mode  can  be  estimated  ignoring  other  contributions.  Rumerman  uses  the  following  relationships. 
When  ,  the  quantity  in  each  of  the  first  two  sets  of  brackets  is  equal  to  unity.  For  Re(v„, ) 

an  integer  n , 

)  =  *■(-  ly  sinh[;r Im(v„  )] « i(- 1)" tt Im(i/„  ) . 

The  wavenumber  pair  isajoint  solution  of  0  —  Z^  Z^(Q.,CX^)—iX^ 

where  XJ!  is  the  reactive  component  of  the  radiation  impedance  and  R'  is  the  resistive  component. 

The  shell  impedance  has  been  taken  to  be  completely  reactive.  At  high  frequencies,  where  the  shell 


36 


impedance  generally  dominates  the  radiation  impedance,  may  be  approximated  by  +iv, , 
where  is  the  real  root  of  0  =  Z*  )-  .  The  imaginary  component  may  be  found  from 

0  =  Zl  {a,  (Q,  aj+i  v,  Zj^  (a,  +  i  v,  Zj^  (Q,  ). 

Therefore,  V,  =  iR'Jzl^  ,  where  =  pcml{K^asmS)  with  cos<5  =  VRl{Ki^a)  and  S  =6^ 
for  the  mode  considered. 

Upon  combining  these  results,  the  pressure  amplitude  backscattered  by  an  isolated  resonant  mode 
into  the  angle  at  which  it  is  traced-matched  with  the  incident  wave  may  be  found.  It  is  simply 


Obviously,  this  expression  is  not  valid  for  arbitrarily  long  shells.  A  limiting  value,  independent  of 
length,  should  be  reached  as  L  — >  <» .  This  limit  can  be  determined  from  the  transformed  expression 
by  noting  that  it  represents  the  scattered  field  as  a  series  of  axial  modes.  When  L  is  very  large,  the 
modal  density  is  very  high,  and  no  one  mode  has  any  significance.  The  solution  must  be  obtained  by 
summing  over  all  modes.  As  L  ,  mTtI L  may  be  thought  of  as  a  continuous  variable  and  the 
summation  approximated  as  an  integral  over  this  variable.  The  integral  itself  can  be  approximated  by 
noting  that  the  integrand  has  singularities  at  values  of  the  integration  variable  at  which  is  an 
integer  n .  The  magnitude  of  the  backscattered  pressure  due  to  the  circumferential  mode  of  order  n , 
and  for  =  Re(ci:„ ),  is  approximately  (Rumerman,  1993). 


""  Rn\\m{a„} 
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MODELING  ISSUES 


APPROXIMATE  SOLUTIONS  OF  THE  DISPERSION  EQUATION 
Formulation 

Recalling  our  results  from  the  section  “Forced  Vibrations  When  p  =  0”  the  system  impedance  in 
that  case  is 


Q.  a 


a^~(Kof 


The  dispersion  equation  =  0  may  be  written  as 


Q  {k,a) 

^  rjr 

2 

pCph  \ 

-Q'  +l-pp.'‘(k,af  +  2(^  +  a\k,aJ 


where  again  k^a=  and  k^=k sin ^  . 

When  p  ^0 ,  the  dispersion  equation  has  a  very  similar  form  (Guo,  1994) 


P\k,ay+i^Z'J\ 

L  ^  _ 

pc^h 

-  O'*  +  — ^  {kfij  +  2(1  +  a'fk^of  -  (1  -  (7^  ]  • 

Note  that  in  the  high-frequency  limit,  the  roots  of  this  equation  are  given  by  the  three  bracketed  terms 
on  the  first  line. 

We  will  follow  Guo’s  approach  (Guo,  1994)  in  taking  these  high-frequency  limits  as  leading  order 
terms,  and  seek  corrections  to  extend  the  solutions  into  the  mid-frequency  region.  Let 


Q' 


Q 


2  ■ 
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The  dispersion  equation  then  takes  the  form 


where  |  is  the  quantity  to  be  determined  and  the  functions  Dq  and  Dj  are-given  by 


1-cr 


^-1 


Ae)=l-^A-i^(|-A)|2CT  +  3-(l-CT=X^-A)] 


with 


fl  = 


a  = 


Po« 

ph 


1 

«  = _ ^ 

^  7 

yOoCo  fel 

It  may  be  seen  that  the  problem  has  been  contrived  in  such  a  way  as  to  suggest  solutions  of  the  form 

Although  this  takes  the  nominal  form  of  an  inverse  power  series  in  frequency,  the  final  results  will 
not  have  such  a  form  because  the  frequency  parameter  also  appears  in  through  both  p  and  ■ 

These  will  be  taken  as  fixed  in  applying  the  expansion.  As  Guo  points  out,  a  strict  inverse  power 
series  in  frequency  is  not  a  suitable  solution  for  cases  of  fluid  loading  because  the  roots  are  generally 
complex  with  real  parts  many  orders  of  magnitude  larger  than  the  imaginary  parts.  An  inverse  power 
series  expansion  would  require  a  large  number  of  terms  in  order  to  obtain  the  first  non-zero 
imaginary  part. 

In  finding  the  roots  corresponding  to  shear  and  compressional  waves,  p  and  are  assumed  to 
be  unaffected  by  the  expansion.  The  reason  for  keeping  p  constant  is  purely  for  convenience.  It  is 
always  much  smaller  than  unity,  under  the  constraint  of  thin-shell  theory,  and  consequently,  makes  a 
negligible  contribution  to  the  roots.  The  reason  for  keeping  z„„  constant  is  subtler.  It  is  because 

Z^  is  the  term  accounting  for  fluid  (radiation)  loading,  and  hence,  the  only  term  that  gives  the  roots 
imaginary  parts.  Keeping  z^n  allow  us  to  find  the  complex  roots  with  only  one 

correction  term.  This  does  not  contradict  the  use  of  the  expansion  above,  because  z^n  of  order 
1/Q ,  as  may  be  seen  from  its  definition  and  our  previous  analysis.  Therefore,  if  an  expansion  is 
begun  at  the  order  1/Q^  ,  z^  should  be  included  in  the  leading-order  expression. 
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For  flexural  waves,  these  arguments  do  not  hold.  Here,  the  dominant  effect  of  fluid  loading  is 
added  mass  with  little  radiation.  contributes  predominantly  to  the  real  part  of  the  roots.  In  fact, 
in  the  frequency  range  of  interest,  the  roots  for  flexural  waves  are  purely  real.  Consequently, 
expansion  of  in  terms  of  frequency  is  required.  Since  flexural  waves  make  a  negligible 
contribution  to  the  scattering  problem,  we  will  not  pursue  their  study. 

To  continue  the  analysis,  we  substitute  the  expansion  into  the  recast  dispersion  equation,  and 
group  terms  of  like  powers  (after  expanding  Dq  in  a  Taylor  series).  To  order  ,  this  leads  to 
solutions  of  the  form 

^  =  A  1  A  (4) 

where  the  prime  denotes  differentiation  with  respect  to  the  argument  and  ^re  the  high-frequency 
limits.  It  may  be  seen  that 

^o(4)=  ^4-1  ^^4-1  + 

and 


Note  that  this  differs  from  Guo’s  expression. 

Compressional  Waves 

The  leading-order  solution  for  compressional  waves  is 

4=1- 

Here, 

^0  (4  )=  1/^  - 1  -  ] 

A(4)=^b+(i-o-)^r 

so  that 
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^_ll  1  [g-  +  (l-cr)Af 

In  terms  of  the  axial  wavenumber,  this  may  be  rewritten  as 


Figure  21  shows  a  comparison  with  the  results  of  Rumerman’s  formulation.  The  parameters  are 
the  same  as  before.  As  can  be  seen,  the  differences  are  minor  and  are  so  for  most  practical  problems. 


0  5  10  t5  20  25 

ka 


Figure  21 .  Comparison  of  approximate  and  “exact”  solution  for 
compressional  modes. 
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Shear  Waves 


The  leading-order  solution  for  shear  waves  is 

2 

1-0 


Here, 


(1-cr)  [4  4 


so  that 


,212  2 

- + 


1  — <7  1/,  \2  1  .A  _V/v^ 

4  4 


This  may  be  rewritten  as 


^  y  2Q^ 

=:; - 

1  — (T 


14(1-..)^  IJi. ' 


1  A  Y 


12^  I  I  1  A  \2  Fi  *  1 

jU--(l-o-)  [1  +  iaz^j 
4 


with 


£2' 


i 

) 

2 

UJ 

-1 

1  . 

•— i- 

2 


2 

(c  ^ 

2 

S 

-1 

Cn 

V  ^  ) 

_ 

K  ^  J 

-l-n 


£2' 


/ 

y 


f 

c 


-|2 


-1 


Note  that  =  0  when 

m 


This  is  the  cutoff  frequency  for  shear  waves.  In  fact,  it  is  the  exact  cutoff  condition  that  can  be 
derived  directly  from  the  dispersion  equation. 

Figure  22  shows  a  comparison  with  the  results  of  Rumerman’s  formulation.  The  differences  are 
insignificant. 
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0  5  10  16  20  25 


ka 

Figure  22.  Comparison  of  approximate  and  “exact”  solution  for  shear 
modes. 


MODAL  WIDTHS 

Previously,  we  had  noted  the  axial  modal  directivity  function 

•  fl  ' 

2  Sin  7;Z 
_  . 
1 

where 


X  =  mn  - kasm^i„^L  . 

As  defined,  B  is  an  even  function  in  x  ^nd  5(o)  =  1 .  Note  that  B  =  ^Iljl  when  A;}f  =  2.004 . 
The  axial  modal  half-width  is  therefore  A(fejtsin  )  =  kxl^  ■  We  will  later  see  that  taking 
l^X  =  1.00  gives  better  agreement  with  observations. 

The  circumferential  modal  directivities  may  be  characterized  in  terms  of  parameters 
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The  circumferential  modal  half- widths  are  then  A(A:flsin^j^  )=  .  Here,  we  do  not  have  an 

explicit  directivity  function.  For  convenience,  we  will  take  the  directivity  to  be  unity  over  the  modal 
width  and  zero  otherwise.  We  will  later  see  that  taking  A^  =  0.25  produces  behavior  in  reasonable 
agreement  with  observations. 

AMPLITUDE 

Our  previous  results  for  the  backscattered  pressure  amplitude,  when  traced-matched  with  the 
incident  wave,  may  be  summarized  as 

This  obviously  ignores  the  transitional  behavior.  )  cannot  be  expressed  in  a  simple  manner. 

However,  in  the  regime  of  interest,  Im(o„  )«=  n  (Rumerman,  1992  and  1993).  The  fai  dependence  is 
less  apparent.  From  our  earlier  results,  we  would  expect 

Ini(or„ )  ~  ■yj {ka  cos  (py  +  {p^a/ phY  j ka  cos (p .  For  typical  problems,  a  good  fit  is 


'  ^  Rn: 


Im(a„  )=  0.366-^- 

Pq  a  ka  cos  (p 

It  is  important  to  note  that,  in  practice,  these  peak  levels  are  not  observed.  There  are  a  number  of 
reasons  for  this.  The  primary  one  has  to  do  with  the  shell  termination.  Caps  modify  the  coupling  of 
the  incident  wave  with  the  shell  modes.  Stiffeners  will  as  well.  Typically,  the  modal  structure  is  not 
significantly  altered,  but  the  backscattered  levels  are  reduced  by  some  coupling  efficiency  factor. 
This  factor  is  usually  in  the  range  of  0.6  to  0.8. 

EXAMPLE 

Let  us  consider  the  example  used  by  Rumerman  (1993)  where  a//i  =  100,  c^/cq  =  3.5  , 
pj Pq  =  7.8 ,  and  o  =  0.3 .  The  model  we  have  presented  produces  the  results  shown  in  figure  23. 

The  color  scale  represents  the  target  strength  in  decibels.  This  may  be  compared  with  Rumerman’ s 
figure  3  (data)  and  figure  7  (predictions).  The  qualitative  agreement  is  very  good.  Given  the 
uncertainties  involved,  the  quantitative  agreement  is  good.  Rumerman  explains  the  reasons  for  many 
of  the  discrepancies.  Figure  24  shows  the  response  over  a  wider  bandwidth. 
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23.  Predicted  response  for  the  example  given  by  Rumerman  (1 993). 


Membrane  Waves 


Aspect  (degrees) 


Figure  24.  Predicted  response  over  a  wider  bandwidth. 
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SUMMARY 


We  examined  the  modeling  of  shell  membrane  waves  and  provided  detailed  theoretical 
background,  both  to  serve  as  a  primer  for  those  unfamiliar  with  the  subject  and  to  set  the  context  of 
the  model  presented.  We  discussed  the  nature  of  the  various  assumptions  and  made  comparisons  of 
the  predicted  behavior  with  “exact”  numerical  solutions,  as  well  as  with  experimental  data.  The 
qualitative  agreement  was  excellent.  Quantitative  agreement  was  good,  given  the  uncertainties 
involved,  as  well  as  sensitivity  to  boundary  conditions. 
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